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Abstract 

We study the effect of a small cutoff e on the velocity of a pulled front in one dimension by 
means of a variational principle. We obtain a lower bound on the speed dependent on the cutoff, 
and for which the two leading order terms correspond to the Brunet Derrida expression. To do so 
we cast a known variational principle for the speed of propagation of fronts in new variables which 
makes it more suitable for applications. 
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I. INTRODUCTION 



In several problems arising in physics, population dynamics, chemistry and other fields, it 
is found that a small perturbation to an unstable state leads to a propagating front joining 
the unstable to an stable state. The simplest model of such phenomenon is provided by the 
scalar reaction diffusion equation 

Ut = U^x + f{u) 

where the reaction term f{u) is a nonlinear function with at least two fixed points, one 
stable and one unstable. Without loss of generality we assume that there is an unstable 
fixed point at m = and a stable fixed point at u = 1. The reaction term f{u) obeys 
additional requirements depending on the phenomenon under study. In the present work 
we shall be interested in two generic classes. The first class, which we label type A, is that 
for which / > in (0, 1), the second class, type B, also called the combustion case, is that 
for which / = in (0, a), and / > in (a, 1). It was proved by Aronson and Weinberger |l] 
that sufficiently localized initial conditions evolve into a monotonic front joining the stable 
to the unstable state. In case B there is a unique speed for which a monotonic front exists. 
In case A, the front propagates with the minimal speed for which monotonic fronts exist. 
This minimal speed satisfies 



2v/f(0) <c* < 2 sup {f{u)/u), (1) 

V 0<M<1 

result also found by Kolmogorov, Petrovsky and Piskunov (KPP) [2]. For the classical 
Fisher-Kolmogorov equation 

Ut = U^x + U{1 - U), 

the upper and lower bounds coincide and the speed is exactly the so called linear or KPP 



value ckpp = 2J f'{0). Fronts for which this is the minimal speed are called pulled since 



this is the speed obtained from linear considerations at the leading edge o 



cases the speed can be calculated from the integral variational principle {4] 

c^ = sup2 



the front. In all 



where the supremum is taken over all positive monotonic decreasing functions g{u) for which 
the integrals exist and where h{u) = ~g'{u). Moreover, the supremum is always attained for 
reaction terms of type B and for reaction terms of type A it is attained whenever c > ckpp- 



Reaction diffusion equations of type A are often used to model phenomena in population 
dynamics, with the assumption that the number of particles or individuals is large. It was 
noticed by Brunet and Derrida ^ that the effect of a finite number of particles can be 
modeled by reaction terms of type A with a cutoff e = where is the average number 
of particles at the saturation state of the front. The effect of a cutoff on the fronts was 
studied for the case f{u) = u — v? and it was found that the speed of the front with a cutoff 
is given approximately by 

(loge)^ 

This result was obtained by a matching approach. Recently this result has been proven 
rigorously by a geometric method [6|]. The precise dependence of the speed of the front on 
the cutoff is not universal. An example of this non universality was shown by introducing 
a small region of vanishing slope next to the cutoff the speed in this case turns out to 
be larger than the KPP value. The effect of a cutoff on reaction diffusion equations is of 
relevance not only as a model of populations with a large but finite number of individuals, 
it is also relevant to the study of noisy fronts [8] and to some problems of particle physics. 
See for example [ol and references therein. 

The purpose of the present work is to show how the speed of a pulled front with a cutoff 
can be found from the variational principle ([2]). It is important to notice that the effect of a 
cutoff of a reaction term of type A is to transform it into a reaction term of type B, reaction 
terms for which the supremum in ([2]) is attained and for which a unique speed exists. 

To obtain this result we reformulate the variational principle in a new way better suited 
to treat the fronts with cutoff. We apply this new form of the variational principle to the 
reaction term considered by Brunet and Derrida, but the results are of more general validity. 
We find a lower bound on the speed which depends on the cutoff, for which the leading order 
is the Brunet-Derrida expression. That is, we show that, 

c > c(e) ^ 2 - + h.o.t 

(loge)^ 

The same result can be obtained by using the alternative variational principle for the speed 

Q. 
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II. A SIMPLER FORM FOR THE VARIATIONAL PRINCIPLE 



In this section we introduce new variables which render the variational formula ([2]) for 
the speed simpler to apply, particularly to the case of fronts with a cutoff. As an application 
of this simpler form we show in the Appendix A how the linear or KPP value ckpp = 2 is 
obtained for the Fisher Kolmogorov equation. 

The variational expression ([2]) implies that for any admissible trial function g{u), 

^2 > 2. ^0 f^^^9{u)du 



lo 9^iu)/h{u)du 

It was shown in \^ that the trial function g{u) for which equality holds diverges at m = 0, 
so it is convenient to consider trial functions which in addition to the requirements g{u) > 0, 
g'{u) < also satisfy g{0) oo. 

Since g{u) is a monotonic decreasing we may perform the change of variables 

u = u{s), where s = 1/g, 

and consider s as the independent variable in ([2]). With this change of variables we find 

/ f[u)g{u)du= +/ — ds, 

Jo So Jo 

where Sq = l/g{u = 1) is an arbitrary parameter and 



Fiu)= r f{q)dq. 
Jo 



du = — d s. 



The denominator becomes 



Iq h{u) Jo \ds 

In this new variables the variational principle becomes 

2 _ ^ Fil)/so + J^'>F{u{s))/s'ds 

C -sup/ 2 , ' W 

u(s) Jo [du/dsj d s 

where the supremum is taken over positive increasing functions u{s) such that u{0) = 0, and 
for which all the integrals in ([3]) are finite. 

In the Appendix A, we illustrate how to use this variational principle to show that for 



profiles satisfying the KPP criterion (i.e., f{u) < f'{0)u, for all < m < 1), c = 2Jf'{0). 
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III. THE SPEED OF THE FRONT WITH A CUTOFF 



In this section we consider the speed of a front for a reaction term with a cutoff. Even 
though we choose a specific reaction term the results obtained are vahd for a larger class 
of reaction terms. We choose the same reaction term studied previously by Brunet and 
Derrida, namely 

\ if < M < e 
u — if e < M < 1. 



For this reaction term 



Fiu) 



if < M < e 

mV2 - m74 - eV2 + eV4 if e < m < 1 



so that F{1) = 1/4 — e^/2 + e^/4. We will show that for a certain trial function the variational 
formula yields 



> c^(e) ^411 



2 
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For the sake of clarity we postpone until section 4 the construction of the trial function. 
Choose the trial function 

{ s if < s < e 

u{s) = \ - - (4) 

I A^/s cos (l){s) if e < s < So 

where 

So = 1/e, (f){s) = ujln{s/e) - (f)^, (5) 

and where 0* is the solution of 

0^,tan0* = 2 1 ^^^1- (''') 

Notice that these definitions imply that in the range e < s < So, —0* < 0(s) < 0*, and that, 
for small e, 0* < 7r/2. Therefore u{s) is positive and monotonic increasing. 

Having chosen a trial function it is straightforward to obtain a lower bound. The numer- 
ator in (I3j) is given by 

so 



N{e) = F(l)/so+ / F{u{s))/s^ds 
Jo 
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6F(1) 



+ — 20, + sin20, - - / —ds 
J 4uj 4 Je 



e2 




+ 




2" 


4y 




e2 




+ 




2 


4 ^ 


4cj 



> eF(l)+ e-- --- +_(20, + sin20,) + O(e3/2(lne) 



(8) 

(See equation ( ]B5I) in the Appendix B for the details on the estimation of the last integral 
in (IH])). The denominator is given by 



Di 



e + 



(du/ds) ds 



We know then that 



l + 4cj')20, + (l-4cu2) sin 20, 
N{e) 



Die) 



(9) 



(10) 



The bound above is rigorous and it is explicitly dependent on e. Expanding c(e) for small e 
we obtain in leading order the desired result, (see part ii of Appendix B, in particular the 
derivation of (IBSP ) 



cHe) =41 



TT 



In el 



+ h.o.t. 



(11) 



It is not difficult to obtain higher order terms in the expansion of c(e) but it is of no 
interest here. 



IV. THE TRIAL FUNCTION 

To find the trial function for which the maximum in the variational formula for the speed 
is atained we should solve the associated Euler-Lagrange equation. This is not possible in 
general since there are few exactly solvable cases. In the present situation we are interested 
in the effect of the cutoff on pulled fronts, that is, on fronts whose speed is determined 
from linearization at the leading edge, therefore we solve the Euler-Lagrange equation in the 
linear approximation. 

The Euler-Lagrange equation for the variational principle ([3]) is 

d'^U ^/('") g 

ds'^ ' 

where A is a Lagrange multiplier. Even though it is unrelated to the present discussion, it is 
worth mentioning that this equation can be obtained by performing the change of variable 
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s = exp(— C2;) in the ordinary differential equation Uzz + cuz + f{u) = 0, hence we identify 
the Lagrange multipher with 

First we obtain the adequate trial function for pulled fronts without a cutoff. In the 
linear approximation the Euler-Lagrange equation is 

subject to m(0) = 0. The solution is of the form u = where a is given by 

a = ^ ± Vl -4A. 

As shown in Appendix A, the best bound is obtained for a 1/2 hence the Lagrange 
multiplier in that limit is A = 1/4. 

Next we construct the appropriate trial function for a pulled front with a cutoff. We 
must solve 

_^ = Mi(0) = for 0<M<e. 

Since the variational formula ([3]) is invariant to scaling in s we may choose, without loss 
of generality, Ui = s. Moreover, since this is valid for < m < e, we conclude that in this 
region < s < e. To sum up we have 

= s, if < s < e. 

For u > e {s > e) but still small enough for the linear regime to be valid, we must solve 

-^ + A— = 0, with U2[e)=ui[e), u^ie) = u^[e) 

The solution to this equation is straightforward, it is given by 

U2{s) = A^/s COS (f){s) 

where 

s . 1 



(j){s) = Lu\n- + 5, tu = -V4A - 1. 

The constants A and S are found matching the solution to ui as indicated above. Applying 
the matching conditions we obtain 

/ 1 _i 
A = ( 1 + ^ j , S = arctan — . (12) 



We must also require that u{s) be positive and monotonic increasing. The first condition 
imphes — 7r/2 < 0(s) < 7i/2. The second condition, that u{s) be monotonic increasing 
imphes 

tan0(s) < — . 

2iu 

Since 0(s) is an increasing function of s we know that 

arctan -— = 0(e) < < (t>{so) < arctan — . 

It is intuitively evident that the best bound will be obtained when the maximum range for 
(p is allowed. We choose then 

0(so) = 0* = arctan (13) 

With this choice, 5 = 0(e) = —0*. The only free parameter left is the arbitrary parameter 
Sq. To fix So we observe that as s ^ Sq the solution must approach u = 1. Since for pulled 
fronts without a cutoff uj = 0, we expect that for small e, u will be small. Then from f|T3|) 
it follows that 

TT 

(l)^ = - - 2uj + h.o.t., 

hence 

m(so) = A^/sqco^^^ ^ A^io sin(2u;) ~ 2ijjA^Js^. 
From f[T2l) we see that for small o;, 

Therefore, for small u;, u(so) ~ y^eso. Requiring that u(so) — > 1 implies then sq = 1/e. With 
this choice for sq it follows that 

0(so) = 0* = —2uj In e — 0,,, 

hence 



a; 



lne| 

Replacing this value of lo in f|T3|) we obtain ([7j), with which the construction of the trial 
function is complete. 
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V. CONCLUSION 



The purpose of this work was to study the effect of a cutoff on the speed of pulled fronts 
making use of the variational formulation for the speed. To do so we have rewritten the 
variational principle in new variables which simplify the problem. An additional advantage 
of this reformulation is that the Euler-Lagrange equation for the maximizer is seen easily to 
be the equation of the traveling front. 

Reaction terms with a cutoff belong to the class of general reaction terms for which a 
maximizer always exists and for which the speed is unique. If the original front without a 
cutoff is a pulled front then, with a cutoff, it is possible to solve the Euler-Lagrange equation 
in the linear approximation and obtain an upper bound on the speed. This value obtained 
from the linear equation is vaid only for sufficiently small cutoffs. The lower bound on the 
speed is a complicated function of the cutoff, the first two terms in the series expansion of 
this bound correspond to the approximate formula found by other approximate means. Here 
we have obtained not only the first two terms in the expansion but a rigorous bound on the 
speed. 

It has been shown that small perturbations of the reaction term close to the cutoff have 
an important effect, and that the Brunet-Derrida term is not universal for all fronts with a 
cutoff. This can be expected since the Euler-Lagrange equation in the linear regime will be 
different in each case. A detailed analysis of this situation will be reported elsewhere. 

In the present work we have studied the effect of a cutoff on a pulled front, the effect of 
a cutoff on pushed fronts and bistable fronts has received less attention. A specific example 



is studied in 



General bounds on the speed have been obtained making use of the 



variationa 
functions 



principle ([2]) and exact solutions have been constructed for piecewise continuous 



12| . These and other related problems will be the subject of future work. 
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APPENDIX A 



Here we show how to recover the value c = 2y/'(0) from the variational formulation ([3]) 
for the speed of propagation of fronts when the profile f{u) satisfies the KPP condition, i.e., 
when f{u) < f'{0)u, for all < u < 1. Let us denote 

F(l) 



N 

So 

and 



+ rFiu{s))\ds, (Al) 

JO 



^0 ( du^"^ 



" = 1 [Ts) "^^ 

For the KPP case, F{u) < /'(0)mV2, hence F(l) < /'(0)/2. Then, it follows from that 

Intregrating the last term by parts and noticing that u{0) = 0, u{so) = 1, and that 
lims^ou/s = u'{0) exists, we get, 

Jo / JO 2 \ as s J zsq Jo s 

Therefore, 

M^<( f'{0)uu'-ds] <f'{OY ^ds {u'{s)fds, (A4) 

VJO S J JO JO 

by Schwarz inequality. However, from (]A3|) we have /'(O) Jq" / (2s^) ds < M, and inserting 
this in flA4p we finally get, 

M < 2/'(0) / {u'{s)yds. (A5) 



Now, from (jA2]), (jAS]), and (jMl) we have that 

2^<4/'(0), (A6) 

for all possible trial functions u. Therefore, taking the supremum of 2N/D over all u, using 
([3]) we finally get. 



c < 2^/'(0). (A7) 
On the other hand, choosing an appropriate maximizing sequence of functions, using the 



variational principle ([3]), we may show that c > 2w/'(0) and thus we can conclude that 



c = 2y/'(0) in the KPP case. For that purpose, just consider the family of trial functions 
Ua = s", which are appropiate trial functions as long as a > |. Evaluating the right side of 
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(j3]) with u = Ua and letting a — > 1/2, we get > 2y /'(O), which combined with (1A7|) yields 
the desired result. Just to illustrate this procedure consider the reaction term f{u) = u — u^. 
Effectively, with this trial function ([3]) implies 

2a - 1 1 2a - 1 2a" 
^ 2 ~ 4(4a- 1)^° _ ■ 

In the limit a ^ 1/2 we obtain > 4. 



APPENDIX B 

In this Appendix we show some details on how to get the lower bound ffTTl) . 

i) Estimating the integral / = / ds for the trial function u given by ([4]). In 

order to estimate this integral we divide it into two parts, as follows, 

h= / -ds, (Bl) 

J e S 

and, 

h= ,^ds. (B2) 
To estimate (IBll) we insert (jlj) and we get 

Ji = cos^ ds < A\t~^/^ - e). (B3) 

As for the second integral we use that m < 1 to get 

/2 < - e. (B4) 

Adding up this two integrals and using ([H]) we see that / can be estimated from above by a 
term of order e^/^| lne|^, i.e., 

/''^'!f!rfs<0(e3/2|lne|^), (B5) 

J e S 

which is small compared with A^/uj e/{Auj'^) = 0(e(lne)'^), since e^/^(lne) — as e — 0. 

We have shown that the contribution of the nonlinear term I can be neglected when 
e — s> 0. To do so we split the integral in two parts. If the reaction term corresponded to 
that of the Fisher- Kolmogorov equation, f{u) = u — v?, a different sphtting is necessary. 
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One can show that the contribution of the nonhnear term always vanishes compared to the 
contribution of the hnear term. A different case is that when the slope of the reaction term 
next to the cutoff vanishes {t], but we are not addressing that problem here. 



ii) Estimating the leading order of ( llOp . From the expressions ([6]) for A and uj we see 
that the leading order in both ([8]) and (Q are the terms proportional to A^/u. Hence, the 
leading order in ( fTTi) is given by 

2(f), + sin 20, + 4uj^ (20, - sin 20, ) ' ^ ' 



which we can write as 



where X = (20, — sin20,)/(20, + sin 20,). Finally, we observe that < X < 1, since 
< 0, < n/2 (in fact, 0, ~ 7r/2), and also that 1/(1 + a) > (1 — a) if a > 0, to conclude 
that 



J > 4(1 - Alu'X) > 4(1 - Alu') ^ 4(1 - - — -). (B8) 

(me)"^ 
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